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Two particle azimuthal correlation functions are presented for charged hadrons produced in 
Au + Au collisions at RHIC (-ysjy„ — 130 GeV). The measurements permit determination of 
elliptic flow without event-by-event estimation of the reaction plane. The extracted elliptic flow val- 
ues {v2) show signiflcant sensitivity to both the collision centrality and the transverse momenta of 
emitted hadrons, suggesting rapid thermalization and relatively strong velocity fields. When scaled 
by the eccentricity of the collision zone e, the scaled elliptic flow shows little or no dependence on 
centrality for charged hadrons with relatively low pr. A breakdown of this e scaling is observed for 
charged hadrons with pr > 1.0 GeV/c for the most central collisions. 

PACS 25.75.Ld 



The primary goal of current relativistic heavy ion re- 
search is the creation and study of nuclear matter at 
high energy densities |l]-||J^. Open questions include 
the detailed properties of such excited matter, as well as 
the existence of a transition to the quark-gluon plasma 
(QGP) phase. Such a phase of deconfined quarks and 
gluons has been predicted to survive for « 3 — 10 fm/c in 
Au -I- Au collisions at the Relativistic Heavy Ion Collider 
(RHIC) 1^, and several experimental probes have been 
proposed for its possible detection and study M. Elliptic 
flow constitutes an important observable ||9|-|l^ because 



it is thought to be driven by pressure built up early in the 
collision, and therefore can reflect conditions existing in 
the first few fm/c. Elliptic flow leads to an anisotropy in 
the azimuthal distribution of emitted particles. A Fourier 
decomposition of this distribution [^7|jl^] 

w^/^d, , « (1 + E 2i^n cos(n(0 - $^))), (0.1) 
d{(l) - <^R) 

provides a characterization of the elliptic flow via the 
second Fourier coefficient V2. Here, is the azimuthal 
angle of an emitted particle and is the azimuth of the 
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reaction plane, defined by the beam direction and the im- 
pact parameter vector . The apparent reaction plane 
is determined from the azimuthal correlations between 
emitted particles, and its dispersion correction from the 
azimuthal correlations between two "subevents" [p^}-pl| . 
An alternative technique for elliptic flow analysis, is the 
Fourier decomposition of the pair-wise distribution in the 
azimuthal angle difference (A0 — (pi — (j)2) between pairs 
of emitted particles [p2-[24|: 



dN 

dA(f) 
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(l + ^2z;2 cos(nA,^)). 

n=l 



(0.2) 



In this case the magnitude of the elliptic flow is char- 
acterize d by the square of the second Fourier coefficient 
in Eq. (0.1), i.e. These methods of analysis can be 
taken as equivalent since (i) the correlation between ev- 
ery particle and the reaction plane induces correlations 
among the particles, and (ii) correlating two subevents 
amounts to summing two-particle correlations The 
results in this letter have several advantages over elliptic 
flow measurements already performed at the same beam 
energy by the STAR [| and PHOBOS § collaborations. 
First, two particle correlations circumvent the need for 
full azimuthal detector acceptance. Second, it allows the 
determination of elliptic flow without event-by-event es- 
timation of the reaction plane and the associated correc- 
tions for its dispersion. Third, this method can serve to 
minimize many important systematic uncertainties (de- 
tector acceptance, efficiency, etc) which may influence 
the accuracy of elliptic flow measurements [ p2|p^ . 

Elliptic flow is predicted and found to be negative for 
beam energies < 4 AGeV and positive for higher beam 
energies in Au -I- Au collisions PJl^-p^ . Recent theoreti- 
cal investigations have made predictions for the centrality 
dependence of the scaled elliptic flow A2 = v^je [ P5|j26| 
where e is the eccentricity or initial spatial anisotropy 
of the "participant" nucleons in the collision zone. The 
eccentricity e shows an essentially linear variation with 
impact parameter 6, for ^flhmax ^ b ^ O.Shmax [D in 
Au -I- Au collisions {bmax ~ 14 fm). For central collisions 
( b^ 5-6 fm), it is predicted that higher energy densi- 
ties are produced and rapid kinetic equilibration in the 
QGP phase leads to a characteristic rise in A2 p5| , p6| . 
In addition, elliptic flow for high px particles has been 
proposed as an observable sensitive to the energy loss of 
scattered partons in a QGP phase [p7| . 

The colliding Au beams (y^s^ — 130 GeV) used in 
these measurements have been provided by the Relativis- 
tic Heavy Ion Collider at Brookhaven National Labora- 
tory (BNL). Charged reaction products were detected in 
the east and west central arms of the PHENIX detector 
P,p8| . Each of these arms subtends 90° in azimuth 0, and 
±0.35 units of pseudo-rapidity rj. The axially symmetric 
magnetic field of PHENIX (0.5 T) allowed for the track- 
ing of particles with pT > 200 MeV/c in the fiducial vol- 



ume of both arms. The Drift Chamber (DC) and a layer 
of Pad Chambers (PCI) located at radii of 2 m and 2.5 m 
respectively, in each arm, served as the primary tracking 
detector for these measurements. A second layer of Pad 
Chambers (PCS), positioned at 5 m in the east arm, was 
employed to confirm the trajectory of charged particles 
which traversed both the DC and PCI. The Zero Degree 
Calorimeters (ZDC), were used in conjunction with the 
Beam-Beam Counters (BBC) to provide the position of 
the vertex along the beam direction as well as a trigger 
for a wide range of centrality selections. 

The present data analysis uses two-particle azimuthal 
correlation functions to measure the distribution of the 
azimuthal angle difference (A0 = (t>i — 4>2) between pairs 
of charged hadrons. Following an approach commonly 
exploited in interferometry studies, a two-particle az- 
imuthal correlation function can be defined as follows 



C(A(/.) 



^cor(A(/>) 

,(A<^)^ 



(0.3) 



where Ncor{A<t>) is the observed A(/) distribution for 
charged particle pairs selected from the same event, and 
Nuncori^ff') is the A(j) distribution for particle pairs se- 
lected from mixed events. Events were selected with a 
collision vertex position, —20 < z < 20 cm, along the 
beam axis. Mixed events were obtained by randomly 
selecting each member of a particle pair from different 
events having similar multiplicity and vertex position. In 
order to suppress an over-efficiency in finding two tracks 
at close angles, hadron pairs within 1 cm of each other 
in the DC were removed from both the Ncor{A(j)) and 
Nuncor{A(f)) distributions. Event centralities were ob- 
tained via a series of cuts in the space of BBC versus 
ZDC analog response |||]; they reflect percentile selec- 
tions of the total interaction cross section of 6.8 barns 
1^. Estimates for the impact-parameter and the eccen- 
tricity, e = (< > — < >)/(< > + < >) were 
also made for each of these selections following the model 
detailed in Ref. ||]. Here, < . . . > represents the spatial 
average (weighted by the density) of participant nucleons 
over the transverse plane of the collision zone 1 1 1 1 . Sys- 
tematic uncertainties associated with the determination 
of e are estimated to be ~ 7%. 

Correlation functions were obtained via two separate 
methods. In the first, charged hadron pairs were formed 
by selecting each particle from a common px range (fixed- 
Pt method). In the second hadron pairs were formed by 
selecting one member from a fixed px range and the other 
from outside this range (assorted-pT method). Within 
statistical errors, both methods of analysis yield similar 
results for the pt range presented below, as would be 
expected for a system dominated by collective motion. 

An important prerequisite for reliable flow extraction 
from PHENIX data is to establish whether or not the 
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180° azimuthal coverage of the detector results in sig- 
nificant distortions to the correlation function. To this 
end, detailed simulations of the detector response, accep- 
tance and efficiency, have been performed for simulated 
data incorporating specific amounts of flow (parameter- 
ized hy V2)- The results from these simulations indicated 
no significant distortion to the correlation functions due 
to the PHENIX acceptance. On the other hand, small 
distortions to the correlation function (for Aip < 25°) as 
well as an incomplete recovery of f 2 could be attributed 
to background contributions. 

These background contributions principally affect the 
extraction of V2 from the correlation function in two ways. 
The distortion to the correlation function at small rela- 
tive angle introduces a small systematic distortion when 
fit with a Fourier function [c.f. Eq. ( |0.2| )]. A good repre- 
sentation of the data was obtained with the fit function 
C{A<P) = A • exp(-O.5(A0/cr)2) + ai • (1 + 2u|cos(2A(/))), 
where the Gaussian term is used to characterize the back- 
ground distortion at small angles. In addition, there is an 
isotropic background of false tracks which are predomi- 
nantly misidentified as high px particles. These contribu- 
tions can be efficiently suppressed in the east central arm 
of PHENIX, by requiring a relatively stringent associa- 
tion between tracks found in the DC and their associated 
hits in PCS. Using such a procedure, the fraction of back- 
ground tracks has been evaluated as a function of pt and 
used to correct V2- Corrections range from ~ 10% at low 
Pt to ~ 25% at 2 GeV/c with a systematic uncertainty 
of 5%. The current analysis is restricted to the range 
0.3 < pt < 2.5 GeV/c to maintain this relatively small 
systematic uncertainty. 

Figures |l|a - d show representative Acj) correlation 
functions obtained for charged hadrons detected in the 
pseudo-rapidity range —0.35 < rj < 0.35. Correlation 
functions for relatively central events (centrality = 20 - 
25%, b ~ 7.0 fm) are shown for hadrons with 0.3 < px < 

2.5 GeV/c and 0.5 < pr < 2.5 GeV/c in Figs. 0a and c 
respectively. The same px selections have been made for 
the correlation functions shown in Figs. 0b and d but for 
more peripheral collisions ( centrality = 40 - 45%, b ^ 

9.6 fm) as indicated. Figs. 0a - d show a clear anisotropic 
pattern which is essentially symmetric about A0 = 90°. 
There is also a visible increase of this anisotropy with in- 
creasing impact parameter and px- These trends are all 
consistent with those expected for in-plane elliptic flow 

©Ojili]. 

The magnitude of elliptic flow and the mechanism for 
its development, can be shown to be related to (a) the 
geometry of the collision zone, (b) the initial baryon and 
energy density developed in this zone, and (c) the de- 
tailed nature of the equation of state for the created nu- 
clear matter |^-|l^,^. Since differential flow measure- 
ments can serve to provide important insights for disen- 
tangling these separate aspects [p^|- p^p5| ], we show the 
results of such measurements in Figs. Wand 0. Fig. S 



shows V2 as a function of centrality for several px selec- 
tions; 0.40 < px < 0.60 (diamonds), 0.60 < px < 1.00 
(squares), and 1.0 < px < 2.5 (circles) GeV/c respec- 
tively. Fig. ^ compares the differential flow V2{px), for 
several centralities as indicated. 

Figs. 1^ and ^ both show relatively large differential flow 
values which increase with increasing impact parameter 
and the px of emitted hadrons. The separate effects of 
spatial asymmetry and the response of the collision zone 
to the generated pressure are also evident in Fig. |3[ That 
is, V2 not only increases with increasing impact parame- 
ter for a flxed px, but also increases with increasing px 
for a fixed centrality selection. Trivially, the magnitude 
of flow should go to zero for very small and very large 
impact parameters. Similarly, its magnitude can be ex- 
pected to be zero for px = 0. It is interesting that Fig. ^ 
indicates an essentially linear rise of V2 with px for each 
of the centrality selections presented. Such a trend can 
not be accounted for via simple geometric considerations 
alone |2|]. However, it is compatible with model calcula- 
tions which assume a strong transverse velocity field [p9| . 
This suggests the presence of strong dynamically driven 
transverse fiow at RHIC. The magnitude and trends for 
V2 shown in Figs. ^ and ^ are consistent with other elliptic 
flow measurements at RHIC H,^. 

Figure aims to disentangle the geometric and dynam- 
ical (px) contributions to the elliptic flow over a broad 
range of centralities or energy densities. To do this, we 
plot A2 as a function of centrality to obtain the dynami- 
cal contributions [^,^ . This evaluation is performed for 
two Px selections (0.40 < px < 0.60 and 1.0 < px < 2.5 
GeV/c) which give rise to relatively low and high px val- 
ues respectively. The underlying idea is that this ratio 
should remove the geometric dependence of V2, while the 
Px selections may provide greater sensitivity to differ- 
ent time scales and energy densities associated with the 
expanding system. 

Figure shows an increase in the magnitude of A2 
with increasing px- This increase can be attributed to 
the dynamical response of the created system, resulting 
from the generated pressure gradients. For hadrons of 
0.4 < Px < 0.6 the observed centrality dependence of A2 
is compatible with s scaling. However, a breakdown of 
this scaling is observed for hadrons with 1.0 < px < 2.5. 
Such a trend may point to a change in the particle 
production mechanism or the possibility that pressures 
larger than those predicted by current hydrodynamic cal- 
culations jl^,|ll| are being produced in the most central 
collisions at RHIC. It is also interesting to note that 
the species composition of the charged particle spectra 
changes dramatically between the two px ranges used in 

Fig.|g. 

To summarize, we have measured two-particle az- 
imuthal correlation functions for charged hadrons pro- 
duced in Au + Au collisions at RHIC (^s„„ = 130 GeV). 
The integral, differential and scaled elliptic flow values 
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extracted from these measurements indicate strong sen- 
sitivity to the colhsion centraUty and the transverse mo- 
menta of emitted hadrons. The centraUty dependence 
of V2 suggests that the high-energy-density nuclear mat- 
ter created at RHIC, efficiently translates the initial spa- 
tial asymmetry into a similar asymmetry in momentum 
space. The pt dependence is consistent with the devel- 
opment of strong transverse velocity fields in the colli- 
sion zone. The centrality dependence of A2 for hadrons 
in the range 0.4 < pT < 0.6 is compatible with e scal- 
ing. However, a breakdown of this scaling is observed for 
hadrons with 1.0 < pr < 2.5. Such a trend could result 
from a number of effects, the most intriguing of which is 
a possible change in the equation of state [g5|j2g]. Ad- 
ditional experimental signatures and model calculations 
will undoubtedly be necessary to test the detailed im- 
plications of these results. Nevertheless, the results pre- 
sented here clearly show that two particle azimuthal cor- 
relation measurements provide an important probe for 
the high-energy-density nuclear matter created at RHIC. 
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FIG. 1. Azimuthal correlation functions for charged 
hadrons as a function of centrality and pr selection. The 
solid curves represent Fourier fits following Eq. (2). Error 
bars are statistical only. 
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FIG. 2. V2 vs. centrality for several pr selections. [F] 
and [A] indicate results obtained with the fixed-pT and 
assorted-pT methods respectively. Systematic errors are esti- 
mated to be ~ 5%; they are dominated by the normalization 
of the correction function for real tracks. For the centrality 
range 0-5% the data points are statistically uncertain and the 
points are omitted. 
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FIG. 3. V2 vs pt for several centrality selections. [F] and 
[A] follow the notation Fig. ^ Systematic errors are estimated 
to be ~ 5%. 
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FIG. 4. The centrality dependence of A2 for two differ- 
ent Pt selections. [F] and [A] follow the notation in Fig. 2. 
Systematic errors are estimated to be ~ 10%, dominated by 
the normalization of the correction function and the model 
determination ofe. 
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